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Abstract 



■ It is proposed that a rearrangement of single-particle degrees of freedom may occur in a 

\ portion of the quantum fluid interior of a neutron star. Such a rearrangement is associated 

with the pronounced softening of the spin-isospin collective mode which, under increasing 
density, leads to pion condensation. Arguments and estimates based on fundamental 
relations of many-body theory show that one realization of this phenomenon could produce 
very rapid cooling of the star via a direct nucleon Urea process displaying a dependence 



2 ' on temperature. 



1. INTRODUCTION 

The EINSTEIN, EXOSAT, and ROSAT orbiting X-ray observatories have measured surface 
temperatures of certain neutron stars and set upper limits on surface temperatures of others 
(0). The data for the supernova remnants in 3C58, the Crab, and RCW103 indicate rela- 
tively slow cooling, while that for Vela, PSR2334+61, PSR0656-I-14, and Geminga point to 
substantially more rapid cooling. In the so-called standard scenario for neutron-star cooling, 
the primary role is played by the modified Urea process {nn — > npe~Ue', npe~ nnue), first 
considered by Bahcall & Wolf (1965) and later reexamined by Friman & Maxwell (1979) in 
terms of an in-vacuum one-pion exchange model. Cooling simulations based on the results of 
these works play out the slow scenario of thermal evolution and fail to explain the rapid cooling 
of some stars. This picture is profoundly altered when in-medium effects are taken into account. 
Modification of A^A^, vrA^, and A'A^ interactions with increasing density may be so strong that 
pion ( [jlT] , p!2t ) and kaon ( prown 1994|| ) condensates form in the interior region of a high-mass 



neutron star, leading to a dramatic increase of the neutrino luminosity (|]T2|, 0, [T^ 
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Here we shall demonstrate that softening of the spin-isospin (pion) collective mode in dense 
neutron matter as predicted by Migdal (1978) could give rise to a rearrangement of single- 
particle degrees of freedom prior to the onset of pion condensation and open a new channel 
of neutrino cooling of neutron stars, by giving access to the direct Urea process. This process 
does not involve a neutron (or other) spectator and, if allowed, is an extremely efficient cooling 
mechanism (||15[|). 

2. BUBBLE REARRANGEMENT OF THE FERMI SPHERE 

A rearrangement of single-particle degrees of freedom takes place if the necessary condition 
for stability of the normal state of a Fermi liquid is violated. At T = this condition requires 
that the change of the ground-state energy Eq remain positive for any admissible variation 
6n{p) of the Landau quasiparticle distribution n{p) away from the normal-state step-function 
distribution 0{p — pF)- Formally, 

5Eo = 1^ ip, n{p)) Hp)-^ > , (1) 

where C,{p,n{p)) = e{p,n{p)) — fi is the energy of a quasiparticle relative to the chemical 
potential /i. The condition (|I]) fails if a depression with ^ < forms in the spectrum ^(p) at 
p > Pf'i it likewise fails if there arises an elevation with ^ > Q aX p < pp- The rearrangement is 
precipitated when the density p reaches a critical value pcF at which there emerges a bifurcation 
and a new root p = po oi the relation 

e(p,n(p);p,^) = 0, (2) 

which ordinarily serves merely to specify the Fermi momentum pp. 

The simplest kind of rearrangement of the momentum distribution n{p) of quasiparticles 
of given spin projection retains the property that its values are restricted to and 1, but 
the Fermi sea becomes doubly connected (0). At densities exceeding the critical value pcF, 
the normal-state distribution 9{p — pp) is modified by the presence of a "bubble," or absence 
of particles, over the range pi < p < pf < pp, with the inner surface pi located relatively 
close to the origin and the Fermi momentum pp readjusted to maintain the prescribed neutron 
density. The distance Pf — Pi between the two new Fermi surfaces lying interior to pp can be 
estimated using the formula ^(p — >■ po, n{p); p) = ^o(p — Pcf) — A{p — poY, where A and are 
positive constants. This formula embodies the essential properties that ^ (p) is negative for any 
p < Pp at p < pcP and that its maximum value first reaches zero at p = pcp. Employing this 
parametrization in the relation (0), it is found that no bifurcation point exists for p < pcP, 
whereas two solutions arise for p > p^p, with the distance between the two new Fermi surfaces 
growing in proportion to y/p — pcp. 

Were such a rearrangement to occur in the neutron subsystem of neutron-star matter, the 
emergence of new Fermi surfaces situated at lower momenta would permit the direct nucleon 
Urea process to operate at a much lower density than hitherto considered possible (|l^, ||). In 
this process, the tandem reactions n — >■ pe~Ue and pe~ nve are driven by thermal excita- 
tions. The condition of high degeneracy prevailing even in young neutron stars implies that 
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these excitations remain close to the Fermi surfaces of the participants. For the direct Urea 
mechanism to contribute appreciably to neutron-star cooling, momentum conservation then de- 
mands satisfaction of the triangle inequalities \pFp — PFe\ < Pfu < Pfp + Pfs among the proton, 
electron, and neutron Fermi momenta. With the conventional singly-connected neutron Fermi 
sphere, these conditions are met only at very high baryon densities where the proton Fermi 
momentum kFp reaches sufficiently large values. The best current estimates yield a threshold 
baryon density of at least five times the saturation density po of symmetrical nuclear matter 
(0). On the other hand, if the neutron subsystem undergoes a rearrangement that allows for 
thermal excitations of neutron quasiparticles at much lower momenta than pFn, specifically at 
Pi and pf in the bubble rearrangement scenario, the triangle inequalities are much more easily 
satisfied and the direct Urea process is greatly facilitated. 

We next provide a quantitative basis for this qualitative idea, by appealing to established 
methods of microscopic many-body theory. For a bifurcation point to arise in the solution of 
equation (H), both the spectrum ^(p) and the scalar component /(pi, P2; k = 0) of the Landau 
amplitude of the quasiparticle interaction must depend strongly on momentum. The connection 
(i, p. 37) 

between these two quantities, wherein ^o(p) is the free spectrum and spin dependence is made 
explicit, is usually employed to find the effective mass M* in terms of the first harmonic of the 
expansion of the amplitude /(pi ■ P2; 0) in Legendre polynomials. We exploit this relation in 
the search for new roots of equation (^. 

Within the integral in (H), the quasiparticle interaction amplitude /(pi,p2;k = 0) may be 
replaced by (M*/M)r*^(pi, P2), the scattering amplitude r'^(pi, P2) = r(pi, P2; k = 0) modified 
by an effective-mass factor ([|ll)- The requisite strong momentum dependences can arise if the 
system approaches a second-order phase transition that occurs at some critical density pc where 
a collective mode of frequency uJs{k) collapses at wave number k = ko and the corresponding 
susceptibility diverges along with the scattering amplitude. Near such a soft-mode critical 
point, the singular part F'^ of the scattering amplitude r(pi,p2;k) at momentum transfer k 
may be expressed quite generally by (0) 

M* 

^K.;i3xiPi^ P2; k, P ^ Pc) = -O^^OfsxDik) + 0„aO/3,D(|Pi - P2 + k|) , (4) 

where O denotes the vertex determining the structure of the collective-mode operator (e.g. 
O = (a ■ n)r for the spin-isospin mode with pion quantum numbers, where n is a unit vector 
along the relevant momentum). In deriving (^), antisymmetry of the two-particle wave function 
under interchange of the coordinates and spins of the two particles has been invoked. The 
propagator D may be parametrized as D{k) = + -^"^{k"^ /k^ — 1)^] where (3{p) measures 
the proximity to the phase transition point, with /3(p pc) ^ Pc — p (cf. Dyugaev (1976)). 

The essential messages of the preceding development are that the singular part r'^(pi, P2; k = 
0) ~ D(|pi — P2I) of the scattering amplitude depends on the difference pi — p2 and that as one 
approaches the soft-mode phase transition point this dependence becomes quite strong. We 
assume that the remaining contributions to r(pi, P2; k = 0) can be adequately incorporated by 
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renormalization of the chemical potential fi. Equation is then easily integrated to produce 
an explicit expression suitable for calculation of the single-particle spectrum, 

m = eip) + 1^0^x0,^) Jd{\p- pii)^(pi)^ • (5) 

We now apply this equation to dense neutron matter in the vicinity of the second-order phase 
transition associated with neutral pion condensation, which is engendered by the softening of 
the spin-isospin mode having 7r° quantum numbers It has been predicted that the 

collapse of this mode will take place at a neutron density pc = Pcw in the range (0.2 - 0.5) fm~'^ 
(roughly 1-3 times po), depending on theoretical assumptions ([|l^, |^). Unfortunately, there is 
as yet no definitive microscopic treatment of neutron-star matter from which one can extract 
or derive quantitatively reliable values for the input parameters f5, 7, and ko of our model. 

In this situation, a reasonable strategy is to perform calculations based on expression (|^) 
for several choices of the parameters of the microscopic model. Substituting (5) into relation 
(2), one finds the critical density pcF for the onset of a bifurcation of the latter equation. For 
P > PcF this equation then determines two new momenta pi and pf where ^{p) vanishes, which 
delimit the bubble region of n{p) and between which ^{p) is positive. Representative numerical 
results for the spectrum ^(p) are plotted in Fig. 1. Results for the phase diagram of dense 
neutron matter in the p/po versus plane are displayed in Fig. 2. Different values of 7 

are considered, while keeping the parameter ko fixed at the value O.QpFn suggested by earlier 
numerical calculations (||12||). 



It is evident that variation of the parameters (3, 7, and ko within sensible bounds can have 
strong effects on the phase diagram and therefore on the extent of the phase with rearranged 
quasiparticle occupation. Nevertheless, our numerical study has documented three salient fea- 
tures of the bubble rearrangement. First, the critical density pcF for the rearrangement is 
less than the critical density pc^ for pion condensation. Since both phenomena stem from the 
strong momentum dependence of the Landau amplitude /(pi,p2;k —>■ 0), rearrangement of 
the quasiparticle distribution may be regarded as a precursor of pion condensation. Second, 
the bifurcation point corresponding to formation of a hole bubble in the neutron momentum 
distribution is positioned at small momenta, po < 0.2pF, irrespective of the applicable value of 
PcTT- Third, the spectrum ^(p) shows a deep depression for p ~ (0.5 — 0.6)pf- And fourth, the 
ratios Pcf/Pot and Po/pf are insensitive to the actual value taken by within the usual range 
of theoretical predictions. 

Analogous considerations apply to the proton subsystem of neutron-star matter, in which 
case one is dealing with the charged pion mode. Estimates ([|TT|, |T^) indicate that this mode 
is also softened in dense matter, with a critical density not far from that for neutral pion 
condensation. One may then argue that, under the influence of the strongly momentum- 
dependent external field provided by the neutron medium, protons will leave the old Fermi 
sphere and occupy states of relatively large momentum, p ~ 0.5pi?„. The impact of this 
further rearrangement on the proton-neutron ratio and on the rate of neutrino cooling requires 
a separate analysis, which we defer. 

Having laid the microscopic basis for a rearrangement of the neutron Fermi surface that 
creates a bubble at low momenta in the Fermi sea, we return to its most striking astrophys- 
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ical implication. Beyond the bifurcation point, the triangle inequalities can now be satisfied 
without the conventional requirement |[)) that the proton fraction exceed some 11-14%. 
Accordingly, the direct Urea process becomes active in the density regime just short of the 
threshold for pion condensation. At temperatures T above the anticipated superfiuid phase 
transition, the T-dependence of the resulting neutrino emissivity is determined through the 
usual expression 

27r 

e{n-^pe~u) = — 20^(1 + 3g\)J26^nn{l - np){l - rie) 

i 

X5(p„ - Pp - Pe - Vv)K^n -Sp-Se- Sy) , (6) 

derived in quasiparticle approximation in terms of the occupations rti and energies Si of the 
reacting particles. We do not consider a renormalization of the coupling constants G"]^ and g\ 
due to medium effects, and the emissivity takes the customary form ([Q) 

e(„ - pe-p) . 1^2 X 10-^^ (^) r» ergcn-sec- . (7) 

where /Xe is the electron chemical potential and the temperature T is measured in multiples of 
10^ K. Since the neutron and proton effective masses M* and M* remain T-independent, the 
bubble rearrangement serves only to turn on the direct Urea process at a lower density, without 
altering the familiar power-law behavior of the emissivity. 



3. FERMION CONDENSATION 



Interestingly enough, there exists a more radical scenario for rearrangement of the quasi- 
particle distribution, known as fermion condensation ([^ ^^)- In this case, the occupancy n{p) 
may be partial, i.e., it may lie between and 1. At T = 0, the new quasiparticle distribution 
n{p) is to be found from the variational condition 

The left-hand side of equation (|[) is just the quasiparticle energy e{p). Hence its coincidence 
with the chemical potential /i means that the spectrum of single-particle excitations will be 
dispersionless (i.e., the quasiparticle group velocity will vanish) throughout the entire momen- 
tum domain VL where ^(p) = e{p) — /i = 0, and not merely at isolated points as implied by 
equation (0). The family of quasiparticles having momenta p G is called the fermion conden- 
sate because of a conspicuous analogy with the low-temperature Bose gas, in which the energy 
of condensate particles is also equal to the chemical potential /x. A key signature of fermion 
condensation has been observed in strongly correlated electron systems (||13|)- flat portions of 
the single-particle spectrum have been seen experimentally in a number of high-temperature 
superconductors ([|17|)- 

The two types of rearrangement - bubble formation and fermion condensation - can com- 
pete with each other in the density regime just below the soft-mode phase transition. Numerical 
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studies ([[T^) demonstrate that fermion condensation wins the contest at nonzero T. Let us sup- 
pose this is the case in the neutron-star medium, while continuing to disregard nucleonic pairing 
phenomena. Consistency of the Fermi-Dirac form n{p,T) = {exp [^(p; n(p, T))/T] + for 
the quasiparticle distribution with the variational condition (^) requires that the spectrum 
^{p, T) of the fermion-condensate phase grows linearly with T at low temperature, implying an 
effective mass inversely proportional to T ([^, Q). 

Although the rearranged momentum distribution derived from equation (|^) differs from the 
bubble configuration, its structure will also admit thermal excitations at low neutron momenta. 
Hence we may again expect most neutron stars to contain a region of relatively moderate density, 
bounded below by pcF and above by pcn, in which the direct Urea process operates vigorously. 
However, due to the new feature of a T-dependent neutron effective mass, M„ cx 1/T, we 
may anticipate an enhancement of the neutrino emissivity relative to the standard result (P]), 
corresponding to a rather than a dependence on the temperature. 



4. CONCLUSIONS 

We have explored the possibility that an effective interaction with strong momentum de- 
pendence gives rise to a rearrangement of the neutron momentum distribution in neutron-star 
matter. Two plausible manifestations of this phase transformation - creation of a doubly- 
connected Fermi surface and fermion condensation - have been considered. Both open the 
prospect that direct nucleon Urea cooling is present in a density regime just below the thresh- 
old for pion condensation and consequently at a density much lower than previously estimated. 
If a fermion condensate is formed, the resulting neutrino emissivity is significantly larger than 
that generated by the direct Urea process in normal matter. Within the affected density range, 
it would therefore dominate all other proposed neutrino cooling mechanisms (|]15|)- Future 



studies along this line will focus on temperatures below the superfluid transition and on the 
effect of the dramatically increased emissivity on neutrino opacity. 
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FIGURE CAPTIONS 



Fig. 1. The dimensionless neutron spectrum |/„(p) = ^nip) / {Pf/'^^) ^-t the critical densities 
PcF corresponding to three different sets of model parameters: (a) 7 = 1.25m7r, fco = ^-^PFm 

= 0.22m2 {p^F ~ 1.19po), (b) 7 = 1.25m^, ko = Q.9pFn, = 0.25 {p,p ~ 1.76 po), (c) 
7 = 1.25m7r, fco = PFn, = 0.13 (pcF — 1-88 po)- Two different positions of the bifurcation 
point, namely po = (for parameter sets (a) and (b)) and po ~ 0.12pi;'„ (for set (c)), are 
indicated by arrows. 

Fig. 2. Phase diagram of neutron matter in the variables p (measured in po) and (mea- 
sured in m^), as calculated for kg — 0.9 ppn and four different values of 7, which (in units) 
label the corresponding the phase boundaries separating the bubble phase (upper left) from 
the normal phase (lower right). 
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